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CHAPTER  1 


Introduction 


The  purpose  of  this  report  is  to  document  the  significant 
results  of  the  first  eight  months  of  our  research  program  under 
contract  No,  DAAD05-76-C-0757.  The  object  of  this  portion  of 
the  program  was  to  gain  an  understanding  of  the  underlying  physics 
of  the  impact  process  such  that  one  might  be  able  to  predict, 
using  a  simple  model,  target  material  damage  under  impact  con¬ 
ditions.  This  knowledge  is  essential  if  one  1b  to  rationally  de¬ 
sign  either  a  more  effective  armor  system  or  a  better  penetrator. 

The  design  of  a  future  penetrator  must  take  into  account  the 
optimum  armor  system  which  may  be  encountered. 

A  useful  tool  for  the  design  of  penetrators  (or  armors)  is 
the  simple  integral  theory  of  impact  which  has  been  developed  at 
A.R.A.P.  during  the  past  several  years.  This  theory  contains  the 
essential  physics  of  the  impact  process,  satisfies  all  of  the 
global  conservation  equations  and  is  contained  in  a  computer  code 
which  is  simple  and  inexpensive  to  operate.  The  integral  theory 
avoids  the  gross  empiricism  of  some  models  and  the  high  cost  and 
complexity  of  multi-element  codes.  Its  simplicity  introduces  a 
degree  of  economy  that  makes  it  reasonable  to  conduct  parametric 
studies.  Availability  of  predicted  trends,  rather  than  single 
point  predictions,  greatly  facilitates  the  interpretation  of 
observations  and  the  selection  of  effective  designs.  The  integral 
theory  can,  therefore,  be  used  to  guide  experimental  programs  and 
to  select  those  designs  which  warrant  the  details  of  the  large 
codes . 

During  the  course  of  our  earlier  studies,  it  was  determined 
that  for  the  purpose  of  calculating  target  response  during  impact, 
it  was  necessary  to  determine  experimentally  at  least  two 
characteristic  quantities  for  any  target  material.  One  of  these 
quantities,  denoted  by  E#  (the  characteristic  energy)  repre¬ 
sented  the  amount  of  energy  required  to  put  the  target  material  In 
a  hydrodynamic  mode.  The  other  quantity,  denoted  by  V#  (a 
characteristic  velocity)  was  a  measure  of  the  elastic  energy  which 
could  be  stored  in  the  target. 

Our  present  results  show  that  an  alternative  pair  of  parameters 
can  be  used  to  shed  more  light  on  the  physics  of  the  impact  process. 
In  this  alternative,  E«  contains  at  least  two  components,  a 
quantity  E#-  which  represents  the  energy  absorbed  during  plastic 
deformation  of  the  target  and  a  quantity  E*e  which  is  the  elastic 
energy  absorbed  by  the  target  during  impact.  E#p  corresponds  to 
the  E#  of  our  earlier  Btudies  and  is  shown  to  be  roughly  constant 
for  a  given  target  material.  E#e  corresponds  to  V#  and  is  shown 
to  be  a  function  of  depth  of  penetration. 

The  work  performed  to  date  and  reported  upon  here  consists 
of  two  parts: 


(1)  experimental  evaluation  of  F.*p  and  F#e  (or 
alternatively  E*  and  V„ )  for  a  broad  spectrum  of  target 
materials,  and 

(2)  theoretical  prediction  of  the  val  "e  or  each  parameter 
using  fundamental  material  properties. 

A  series  of  impact  tests  was  conducted  in  the  A.R.A.P. 

Impact  Facility  on  sixteen  target  materials.  The  penetration 
data  obtained  from  these  tests  were  used  in  conjunction  with  the 
simple  code  of  the  integral  theory  to  evaluate  the  characteristic 
properties  for  each  target  material.  Simultaneously,  a  theory 
has  been  developed  which  relates  the  characteristic  properties 
to  fundamental  material  properties,  such  as  hardness  and  elastic 
modulus,  which  can  be  measured  in  static  tests.  This  theory 
accounts  for  both  strain-rate  effects  and  material  property  change 
due  to  shear  heating  in  the  deformation  region.  Because  it 
clearly  identifies  the  fundamental  properties  which  most  in¬ 
fluence  material  response,  the  theory  can  be  used  to  identify 
candidate  armor  (or  penetrator)  materials. 

In  what  follows,  we  will  briefly  review  the  integral  theory 
of  impact  in  Chapter  2.  The  experimental  data  which  have  been 
analyzed  to  date  and  the  data-theory  correlations  will  be  shown 
in  Chapter  3.  The  development  of  the  theory  for  E#p  and  E*e 
will  be  given  in  Chapter  4.  Finally,  a  summary  of  conclusions 
will  be  given  in  Chapter  5. 


CHAPTER  2 


Integral  Theory  for  a  Rigid  Sphere 

In  this  chapter  we  will  outline  the  equations  which  are  used 
to  evaluate  E-  and  V*  for  a  typical  target  material.  The 
integral  theory  of  impact  starts  with  three  basic  equations  which 
govern  the  partition  of  energy  between  the  center  of  mass  kJnetic 
energy  of  an  impacting  projectile  and  the  work  done  on  the  target 
and  on  the  projectile  during  impact.  These  equations  take  a  vejy 
simple  form  for  the  normal  impact  of  an  incompressible  cubic 
pellet  that  remains  a  rectangular  parallelipiped  during  impact. 
The  equations  defining  the  motion  and  deformation  of  such  an 
idealized  projectile  have  been  described  in  earlier  reports 
(Refs.  1  and  2)  and  will  not  be  repeated  here.  Instead,  we  will 
concentrate  on  the  penetration  of  a  nondeforming;  sphere  -  the 
geometry  which  was  studied  in  the  test  program. 

Figure  1  schematically  depicts  the  penetration  of  a  semi¬ 
infinite  target  by  a  nondeforming  sphere.  The  momentum  equation 
for  the  sphere  may  be  written 


§t(mpV  ■  - p 


(1) 


where  F  is  the  total  drag  force  imposed  by  the  target  on  the 
projectile,  mD  is  the  projectile  mass,  and  Vc  is  the  center  of 
mass  velocity.  In  this  report,  we  will  consider  only  normal  im¬ 
pacts.  Hence,  the  velocity  vector,  Vc  ,  will  always  be  normal 
to  the  surface.  Because  the  projectile  is  nondeforming  and  is  not 


V  at  every  point  of  the 


spinning,  the  axial  velocity  is 
projectile . 

If  both  sides  of  (1)  are  multiplied  by  Vc 
an  equation  for  the  particle  kinetic  energy, 
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the  result  is 


(2) 


The  rate  at  which  work  is  done  on  the  target  is  given  by: 

3twt-Fvc  (3’ 

The  addition  of  (2)  and  (3)  shows  that  for  a  nondeforming  pro¬ 
jectile,  the  initial  kinetic  energy  can  only  be  converted  into 
work  done  on  the  target. 

The  energy  transferred  to  the  target  is  distributed  between 
kinetic  energy  and  other  nonkinetic  forms  of  energy  depending  on 
the  state  of  the  target  material.  This  can  be  written: 
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For  the  rigid  sphere,  (4)  may  be  written: 
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where  Is  the  target  mass  density,  VA  Is  the  local  velocity 

normal  to  the  surface  of  the  sphere  and  Is  a  function  of  the 
angle  $  from  the  nose,  and  A  Is  the  submerged  or  Interface 
area  which  Is  also  a  function  of  $  .  The  integration  is  taken 
over  the  entire  submerged  area  of  the  sphere.  This  equation  may 
be  taken  as  a  definition  for  both  Cd  and  E#  and  we  may 
roughly  identify  Cp  with  the  hydrodynamic  drag  of  the  target 
and  E„  with  the  energy  required  to  put  the  target  naterial  into 
a  hydrodynamic  mode. 


If  we  equate  (3)  and  (5),  we  obtain  an  expressicn  for  the 
total  force  at  the  interface: 
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We  see  that  the  force  consists  of  two  components.  The  first 
component  is  caused  by  the  acceleration  of  target  material  around 
the  sphere.  This  term  dominates  when  Vj  >>  E#  ;  i.e.,  target 
inertia  slows  the  projectile  when  the  velocity  is  large.  The 
second  component  is  independent  of  location  on  the  sphere  and  is 
due  to  energy  storage  in  the  target  material.  This  term  dominate 
when  £  H  >>  , 


If  we  substitute  (6)  back  into  (1),  we  get  the  i_quation  of 
motion  for  the  sphere: 
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we  get  the  differential  equation: 


where 

depth, 


2  1 
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<t>s  ,  the  submergence  angle,  is  related  to  the  penetration 
p  ,  by  the  equation: 


cos  $ 


1  “  P/R 

If  we  perform  the  area  integration  in  (11),  we  obtain: 
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Substitution  of  (12)  into  (13!  yields  a  relation  between  Vc  and 
p  ,  A  second  equation  between  these  two  parameters  is  given  by 
the  definition  of  velocity: 


jL 

dt 
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(14) 


Equations  (12)  through  (14)  can  be  integrated  simultaneously 
to  yield  p(t)  and  Vc(tj  .  The  integration  proceeds  from  , 
t  *  0  to  t  ■  t„  ,  the  time  at  which  Vc  is  reduced  to  V„  . 

Three  parameters  appear  in  these  integral  relations ;  CD  , 

E,  and  V#  .  If  one  is  known,  then  the  equations  can  be  used  to 
evaluate  the  other  two  for  a  given  impact  velocity  and  crater 
depth.  The  hydrodynamic  drag  coefficient  for  a  sphere  is  of  order 
unity  and  varies  within  relatively  narrow  bounds.  In  our  previous 
studies  of  deforming  particles  in  Dupersonie  flowfields  (e.g.,  water 
drops  in  shock  layers)  we  have  obtained  good  correlation  between 
theory  and  data  using  Cp  ■  2  .  This  value  can  be  obtained 
using  the  Newtonian  approximation  for  the  pressure  induced  on  the 
surface  of  a  grossly  deformed  projectile  -  the  limiting  shape 
being  a  flat  disk.  Newtonian  theory,  to  first  approximation, 
states  that  the  force  induced  on  the  surface  is  due  to  the  de¬ 
struction  of  the  normal  component  of  momentum.  For  a  flat  disk, 
the  normal  momentum  is  <;  and>  therefore, 


FDRAG 


ptVcAdisk 


(15) 


If  we  substitute  (15)  Into  the  equation  which  usually  defines  the 
drag  coefficient  _  CD  -  FCRAG/ptcA/2  we  See  that  CD.  “  2  1 

If  we  apply  the  Newtonian  approximation  to  a  rigid  sphere, 
we  find  that  the  drag  force  is  given  by: 


^An  alternative  approach  in  which  the  integration  proceeds  to 
Vc  *  0  and  an  elastic  erergy  term  is  included  in  (13)  is 
discussed  in  Chapter  4. 
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The  Newtonian  approximation  is  valid  only  at  the  outer  edge 
of  the  disturbed  region  and  not  at  the  projectile  surface, 
although  this  distinction  is  generally  ignored.  To  obtain  the 
pressure  at  the  surface,  a  centrifugal  force  term  must  be  in¬ 
cluded  in  the  momentum  equation  to  account  for  curvature  effects. 
This  term  has  the  net  effect  of  reducing  CD  to  approximately 
0. 7b  for  a  sphere .  u 

The  analysis  which  is  discussed  in  the  next  chapter  web 
based  on  Cp  *  2  .  In  the  future,  the  analysis  will  be  repeated 
using  the  more  appropriate  value  of  Cp  for  the  sphere. 

In  either  case,  with  Cp  ■  constant,  we  can  return  to  (13) 
and  (1*0  and  numerically  solve  these  equations  to  obtain  the 
combination  of  V#  and  E#  which  best  matches  the  data  for  each 
target  material. 


CHAPTER  3 


Experimental  Results 

This  chapter  describes  the  A.R.A.P.  Impact  Facility, 
summarizes  the  test  program,  presents  the  experimental  data  and 
the  evaluation  of  E*  and  V„  for  each  of  the  target  materials. 

A.R.A.P.  Impact  Facility 

Figure  2  depicts  schematically  the  A.R.A.P.  Impact  Facility. 
This  facility  consists  of  a  mounted  weapon,  and  enclosed  test 
tube  and  test  chamber.  The  weapon  used  for  most  of  the  tests 
is  a  Winchester  270  caliber,  smooth  bore  rifle  which  is  perman¬ 
ently  mounted  to  a  support  and  is  bure-sighted  on  the  target. 
Cartridges  are  hand-loaded  using  Hercules  2400  gunpowder  and  the 
rifle  is  remotely  fired. 

The  projectiles  are  0.250-inch  diameter  balls.  All  testing 
to  date,  using  0.250-inch  diameter  balls,  has  been  done  with 
tungoten  carbide  balls.  Future  tests  will  include  several  other 
ball  materials  including  lead,  aluminum,  and  glass.  The  balls 
are  mounted  at  the  end  of  the  cartridge  using  a  bore-fitting 
Lexan  sabot.  For  tungsten  carbide  balls,  the  velocity  range  of 
the  rifle  is  700  to  5,000  feet  per  second. 

The  sabot  is  manufactured  in  sections.  It  separates  aero- 
dynamically  upon  leaving  the  nozzle  and  each  piece  hits  the 
stripper  plate  located  at  the  downstream  end  of  the  test  tube. 
Only  the  projectile  enters  tne  test  chamber  and  hits  the  target. 

The  target  is  mounted  in  a  permanent  holder  attached  to  the 
downstream  end  of  the  test  chamber.  Most  of  the  targets  have 
been  circular  disks  with  a  nominal  diameter  of  6  inches  and  a 
thickness  of  1  inch.  However,  the  holder  is  versatile  and  can 
accommodate  any  shape  which  can  fit  within  a  6-lnch  diameter 
circle . 

The  velocity  of  the  ball  is  measured  UBing  a  Schmidt-Weston 
Chronograph.  Two  light  screens,  two  feet  apart,  sense  the  pass¬ 
age  of  the  ball  using  a  photo  resistor  element.  The  flight  path 
of  the  ball  is  illuminated  by  light  which  passes  through  slits 
in  a  shelf  in  the  test  chamber.  The  shadow  produced  by  the  ball 
on  the  first  screen  triggers  a  counter  and  a  shadow  on  the 
second  screen  stops  the  counter.  A  digital  readout  of  the  velo¬ 
city  is  provided  on  the  display  board. 

In  order  to  extend  the  low  velocity  range  of  the  facility, 
we  also  use  a  Power  Line  880  Air  Oun.  This  rifle  has  a  velocity 
range  between  160  and  740  feet  per  second  when  firing  0.156  and 
0.172  inch  diameter  tungsten  carbide  and  chrome  steel  balls. 

The  lower  limit  of  velocity  is  set  by  the  chronograph. 
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For  each  target  material,  tests  were  also  conducted  by 
dropping  a  ball  from  the  laboratory  ceiling  onto  the  target. 

From  a  height  of  approximately  14  feet,  the  impact  velocity  was 
computed  to  be  28  feet  per  second. 

Test  Program 

Table  1  summarises  the  status  of  the  test  program.  A 
total  of  141  tests  have  been  conducted  on  16  target  materials. 
These  targets  include  7  nearly  pure  metals:  aluminum,  cadmium, 
copper,  lead,  iron,  silicon,  and  zinc;  3  alloys;  mild  steel, 
hard  armor  steel  and  armor  aluminum;  and  a  range  of  ceramics 
and  composites,  including  glass,  acrylic,  polycarbonate,  sodium 
chloride,  boron  carbide  and  Kevlar. 

Table  1  shows  the  number  of  tests  and  the  velocity  range 
for  each  target  material.  The  lower  limit  of  the  velocity 
range  is  near  the  chronograph  limit.  The  upper  limit  is  set  by 
a  criterion  which  limits  the  crater  depth  to  half  the  target 
thickness  in  order  to  minimize  backface  effects. 

The  status  of  the  test  program  is  shown  by  the  three  columns 
on  Table  1.  The  first  column  shows  that  all  of  the  test  targets 
have  been  fabricated  except  for  silicon  and  boron  carbide. 

These  targets  will  be  completed  shortly.  Testing  has  been 
completed  for  each  of  the  target  materials  except  the  above- 
mentioned  two  and  Kevlar. 

The  target  material  properties  E*  and  V#  have  been 
evaluated  for  the  nine  materials  which  are  checked  in  the  analysis 
column  of  Table  1.  These  results  are  described  below.  The 
remaining  materials  will  be  completed  in  the  near  future. 

Data  Analysis 

Figure  3  shows  the  penetration  data  for  a  soft  aluminum 
(1100-F)  target.  The  figure  presents  maximum  crater  depth 
normalized  by  projectile  diameter  versus  projectile  Impact  velo¬ 
city.  Data  are  shown  for  both  tungsten  carbide  and  steel  ballB. 
The  impact  velocity  varied  from  28  feet  per  second  to  1700  feet 
per  second.  For  the  latter  velocity,  the  ball  remained  embedded 
in  the  target  and  the  maximum  crater  depth  was  slightly  in 
excess  of  half  the  target  thickness.  In  all  of  the  tests,  there 
was  no  measurable  plastic  deformation  of  the  ball  after  impact. 

The  solid  lines  show  the  computed  penetration  using  the 
equations  discussed  in  Chapter  2  and  the  model  parameters 
E„  ■  79  BTU/lbm  and  V„  ■  0  ft/sec.  The  correlation  between 
theory  and  data  is  very  good  for  a  two-decade  range  of  impact 
velocity  and  a  factor  of  two  range  in  density  ratio.  All  of 
these  computations  are  based  on  Cp  «  2  .  The  effect  of  using 
a  more  appropriate  drag  coefficient  for  a  rigid  sphere  (Cp  ■  I) 
will  be  described  at  the  end  of  this  chapter. 
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Figure  4  shows  the  experimental  value  of  Eg  for  each 
test  condition  as  a  function  of  the  normalized  penetration. 

This  value  of  Eg  is  obtained  numerically  by  matching  the 
measured  penetration  with  theoretical  computations  for  each 
impact  test.  It  can  be  seen  that  the  data  fall  within  a  narrow 
band  and  that  Eg  is  approximately  constant.  The  one  excep¬ 
tion  to  this  result  is  a  shallow  penetration  data  point  which  is 
relatively  Inaccurate.  The  solid  line  shows  the  average  value 
of  Eg  ■  79  BTU/lbm  and  it  is  this  value  which  was  used  in  the 
theoretical  computations  in  figure  3. 

The  dashed  curve  shows  the  present  theoretical  prediction 
for  Eg  .  This  theory  is  described  in  the  next  chapter.  For  the 
present,  note  that  the  upward  slope  of  the  curve  for  decreasing 
p/do  is  the  elastic  contribution  to  Eg  .  The  asymptotic  value 
for  large  p/do  is  the  plastic  contribution  to  Eg  . 

Figure  5  shows  the  penetration  data  for  pure,  open  oast 
cadmium.  The  correlation  between  theory  and  data  is  good  for 
velocities  below  1,000  feet  per  second  using  Eg  ■  27  BTU/lbm 
and  Vg  *  0  ft/sec.  For  larger  velocities,  the  theory  under- 
prediots  the  crater  depth.  For  these  velocities,  the  dominant 
contribution  to  the  drag  force  is  inertia  of  the  target  material. 
Beoause  we  used  CD  ■  2  rather  than  cn  *  “1  for  these  computa¬ 
tions,  we  have  overpredicted  the  drag,  and,  therefore,  under- 
predioted  the  penetration.  For  lower  velocities,  Eg  is  the 
dominant  factor  and  Cn  has  only  a  slight  effect.  Examples  of 
the  effect  of  CD  are  given  at  the  end  of  this  chapter. 

Figure  6  shows  the  computed  values  of  Eg  .  Again,  because 
of  the  CD  effect,  the  theory  overpredicts  the  value  of  Eg 
for  large  penetrations.  It  will  be  shown  later  that  a  lower  value 
of  CD  will  tend  to  increase  the  computed  value  of  F*  for  the 
large  p/do  .  To  a  first  approximation,  the  average  value  of  Eg 
for  cadmium  is  27  BTU/lbm. 

Figures  7  and  8  show  the  data  and  theoretical  computations 
for  hot  rolled  electrolytic  tough  pitch  copper.  The  correlation 
between  theory  and  data  is  good.  To  a  first  approximation 
Eg  ■  35  BTU/lbm  and  V,  «  0  ft/sec. 

Figures  9  and  10  show  the  data  and  theoretical  computations 
for  Class  40  gray  cast  iron.  Good  correlation  between  theory 
and  data  is  obtained  using  Eg  ■  133  BTU/lbm  and  Vg  ■  25  ft/sec. 

« 

The  results  for  pure,  open  cast  lead  targets  are  shown  in 
figures  11  and  12.  Lead  has  the  lowest  value  of  Eg  of  all  the 
materials  tested.  Hence,  the  drag  force  on  the  projectile  is 
primarily  due  to  inertia  of  the  lead  target.  We  have  already 
concluded  that  CD  ■  2  is  inappropriate  for  a  rigid  sphere  and, 
therefore,  it  is  not  surprising  that  the  theory  underpredicts 
the  penetration  for  large  velocities,  the  regime  where  inertial 
effects  are  dominant.  Note  that  some  of  the  lead  targets  were 
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2.5  inches  thick.  Hence*  open  symbols  are  shown  for  deeper 
penetrations  than  the  closed  symbol  in  figure  11. 

Figures  13  and  1^  show  the  results  for  mild  steel  (H.R.  1020), 
and  figures  15  and  16  are  for  pure,  open-cast  zinc.  The  correla¬ 
tion  between  theory  and  data  is  good  for  each  metal.  For  steel, 

E#  *  165  BTU/lbm  and  V*  ■  50  ft/sec.  Note  that  the  rigid 
particle  correlation  is  also  good  for  the  projectiles  which 
fractured  upon  impact.  This  is  because  the  tungsten  carbide 
ball  has  had  very  little  deformation  prior  to  brittle  fracture. 

For  zinc,  E#  ■  58  BTU/lbm  and  V#  s  o  ft/sec. 

The  final  material  for  which  data  analysis  is  complete  1b 
sodium  chloride.  Salt  is  the  most  brittle  of  the  targets 
discussed  in  this  chapter.  Figure  17  shows  very  good  correlation 
between  theory  and  data  for  E*  ■  82  BTU/lbm  and  V#  *  0  ft/sec. 

The  computed  values  for  E#  are  shown  in  figure  18.  There  is 
more  scatter  in  the  data  which  is  to  be  expected  from  the  brittle 
nature  of  the  target  material.  However,  most  of  the  data  do  lie 
in  a  band  between  E#  ■  70  BTU/lbm  and  E*  ■  90  BTU/lbm. 

Table  2  summarizes  the  target  material  properties.  The 
experimental  and  theoretical  values  of  E#  are  contained  in  the 
last  two  columns.  The  basis  for  the  theoretical  values  is  de¬ 
scribed  in  the  next  chapter.  The  experimentally  obtained  values 
of  Em  are  shown  for  the  8  materials  described  above,  based  on 
C_  ■  2  and  for  two  materials  (polycarbonate  and  silicon)  for 
which  data  analysis  is  not  complete.  The  values  in  parentheses 
are  the  E#  values  for  copper  and  lead  based  on  CD  ■  1  . 

The  agreement  between  theory  and  data  for  E#  is  excellent. 
Note  that  preliminary  results  for  typical  armor  materials  show 
that  E*  *  200  BTU/lbm  .  Also,  note  that  a  material  has  been 
tested,  boron  carbide,  which  exhibits  E#  *  1200  BTU/lbm  .  The 
basis  for  these  calculations  and  their  implication  is  described 
in  the  next  chapter. 

Cp  Effects 


We  have  noted  that  all  of  the  computations  described  above 
were  based  on  a  constant  value  of  the  hydrodynamic  drag  coef¬ 
ficient  appropriate  for  a  deformed  projectile  (Cp  ■  2)  .  A  more 
appropriate  value  for  a  rigid  sphere  is  CD  ■  1  .  We  have 
investigated  the  effect  of  CD  for  two  materials.  Figure  19 
shows  the  data-theory  correlation  for  the  copper  target.  The 
dashed  curve  is  for  C_  »  2  and  is  reproduced  from  figure  7. 

The  solid  curve  is  foru  CD  ■  1  and  results  in  much  better 
correlation  with  the  datat  To  first  approximation,  the  value  of 
E#  increases  about  9X  when  Cn  is  reduced  from  2  to  1. 


Figure  20  shows  the  computed  value  for  E#  for  Cp  ■  1. 
This  curve  (compared  with  figure  8)  shows  that  E„  is  nearly 
constant  for  large  p/do  . 
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COPPER  TARGET 


Figures  21  and  22  show  the  effect  of  Cp  for  lead.  Again, 
much  better  correlation  with  data  is  obtained  for  Cp  ■  1  .* 

In  this  case,  the  value  of  E„  is  increased  by  nearly  5055  with 
the  lower  value  of  CD  . 


4*  '  ■  -  ,  ------ 

The  data  suggest  that  Cp  may  be  slightly  less  than  1  .  There 
is  theoretical  justification  for  a  lower  value  of  drag  coefficient. 
Indeed,  when  centrifugal  force  effects  are  included  in  the 
Newtonian  pressure  approximation  Cp  *0,75  . 
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CHAPTER  4 


Theory  for  Characteristic  Properties 

The  total  initial  kinetic  energy  of  a  nondeforming  pene¬ 
trator  is  absorbed  by  the  target  in  three  different  forms 
(Eq,  (16)).  Part  of  the  energy,  given  by  goes  into  accel¬ 

erating  the  target  material  as  it  moves  around  the  penetrator. 

It  is  accounted  for  by  the  drag  term  CpV2  in  the  deceleration 
equation.  This  term  has  been  discussed  above.  The  second  term, 
Up  ,  is  the  energy  that  goes  into  plastic  deformation  of  the 
target  which  is  rapidly  dissipated  as  heat.  The  third  term, 
ue  ,  accounts  for  the  elastic  or  recoverable  energy  absorbed  by 
the  target  in  elastic  deformation.  Thus,  we  may  write 

f  Vo  •  KDRA0  *  up  *  Ue  <16> 


Each  of  these  terms  can  be  predicted  in  advance  from  the 
hardness,  Young's  modulus,  melting  temperature  and  other  quantities 
obtainable  from  static  tests.  In  this  section  a  theoretical 
analysis  is  piesented  leading  to  formulas  for  Ue  and  Up  . 

These  formulas  are  then  used  to  predict  in  advance  the  behavior 
of  armor  aluminum,  armor  steel,  and  boron  carbide,  First,  we 
shall  derive  Ue  and  then  derive  Up  . 

Theory  for  Elastic  Energy 

There  are  two  approaches  to  including  the  elastic  energy 
in  the  integral  formulation  of  impact: 


(1)  The  dissipative  work  done  on  the  target  due  to  drag 
and  plaatic  deformation  can  be  calculated  up  to  the  point  where 
the  penetrator  velocity  Vp  <  V#  ,  where  V„  represents  the  ve¬ 
locity  of  the  particle  at  which  all  its  remaining  kinetic  energy 
can  be  absorbed  elastically.  Then 


V. 


DEPTH  OF  PENETRATION  p 


■/  Vdt  ‘ '/ 

V.  V. 


Alf( 


Vm  dV 

■  J?.  ■ 


(17) 


2 


+  E 


The  elastic  energy  then  appears  as  an  integration  cutoff  or  an 
effective  constant  of  integration. 


(2)  Alternatively,  one  can  treat  the  elastic  energy  as  a 
volume  work  of  the  same  type  as  Ep  .  Then  the  elastic  energy 
can  be  brought  inside  the  deceleration  equation: 


in  V  »  mV 
P  P 


d V 
dp 


(18) 
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Since  the  elastic  energy  is  a  surface  effect,  we  expect  it 
to  be  a  function  of  p/d  .  The  velocity  cutoff  when  integrating 
penetration  in  this  case  will  be  zero: 


To  calculate  the  elastic  energy,  assume  some  plastic  work  has 
already  been  done  to  produce  a  depression  in  the  target  of  depth  p 
(see  fig.  23a).  For  a  nondeforming  ball  seated  in  such  a  de¬ 
pression  we  may  assume  the  elastic  behavior  of  the  target  is 
linear.  If  r  is  the  radius  of  the  contact  area,  d  the  ball 
diameter,  then  a  force  F  on  the  ball  produces  an  elastic  stress 
P/irrs  and  an  elastic  strain  k^Cx/r),  where  ki  is  a  constant  re¬ 
flecting  the  average  strain  over  the  deformation  field  (fig.  23b): 


PRESSURE 


■i(S) 


E  is  Young's  modulus.  k1  shall  be  obtained  from  static  tests. 

The  maximum  elastic  force  F#  the  target  can  sustain  without 
plastic  deformation  is  given  by  its  Brinell  hardness: 

F„  ■  irr2B  (20) 

2 

assuming  the  contact  area  irr  doesn't  change  during  loading. 


Then  the  total  elastic  energy  which  can  be  absorbed  is: 

x#f#  mV# 


-F 


x#  can  be  obtained  from  (19),  yielding: 

F* 

x»  "  iTrEFT  (2t 

and  1 

Ue  -  B2r3ir/2Ek1  (2; 

Using  the  relation  between  contact  radius  r  and  depth  p 
for  a  spherical  indenter  of  diameter  d  , 


/cip  - 


as  well  as  (21), one  may  obtain  an  expression  for  V# 


I  p  3 

bVV 

m  Ek. 

P  1 


B2d3it 


f  (D'Vl)'” 


37 


Figure  23 


where  ki  is  to  be  found  from  static  tests,  y  is  the  strain 
rate  factor  to  account  for  the  increase  in  Brinell  hardness  at 
high  strain  rates.  The  strain  rate  hardening  of  B  may  be 
obtained  from  the  velocity-temperature  relationship: 


B(T, e )  -  B(T( 1  -  61n(e/e0)) 


(26) 


which  equates  high  strain  rates  with  lower  temperature,  combined 
with  the  temperature  dependence  of  B,.  Typically,  Brinell  hard¬ 
ness  is  measured  at  strain  rates  of  c  -  10"‘/sec  ,  while  impact 
tests  have  e  *  10*/sec  .  This  leads  to  typical  values  for  y  of 
1.5  for  metals  and  5-10  for  some  plastics. 

From  the  point  of  view  of  interpreting  impact  data,  we  have 
found  it  more  convenient  to  convert  the  expression  for  elastic 
energy  into  an  energy  per  volume,  or  energy  per  mass,  E«e  ,  where 
the  volume  taken  is  that  of  the  indentation.  This  is  the  second 
approach  to  treating  the  elastic  energy.  It  is  then  directly 
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where  the  denominator  is  the  volume  of  a  spherical  section.  This 


simplifies  to: 
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In  order  to  obtain  the  parameter  kj  of  the  elastic  de¬ 
formation  field  and  to  check  the  predicted  (p/d)  dependence  of 
V„  and  E„e  empirically,  we  have  carried  out  a  series  of  static 
tests.  Using  a  Hounsfield  tensometer,  which  is  capable  of 
applying  up  to  a  2-ton  compressive  load,  a  tungsten  carbide  ball 
was  pressed  into  several  target  materials  and  the  elastic  energy 
Ue  measured  by  the  following  technique: 

(1)  Indent  target  to  depth  p  using  tungsten  carbide  ball 
and  Hounsfield  tensometer  at  some  load  G  . 

(2)  Back  out  ball  and  measure  depth  p  . 

(3)  Reinsert  ball,  reapply  compressional  load  slowly, 
measuring  distance  of  travel  on  vernier  micrometer. 

(l|)  Periodically  back  out  ball  and  check  depth  to  make  sure 
it  hasn't  increased. 


(5)  Measure  maximum  force  which  can  be  applied  without 
producing  additional  plastic  deformation.  This  is  F#  . 

(6)  The  corresponding  length  of  travel  of  ball  is  x#  . 

(7)  u  .  ZjXjl 


Using  the  value  of  Ue  obtained  by  this  technique*  we  may 
obtain  values  for  V*  and  E*e  from  the  static  tests  using  the 
expressions! 

V*  ■ 


E#e  "  p^TOTT  (29) 

where  y  is  the  strain  rate  factor  introduced  above. 

The  empirical  values  of  V#  obtained  in  this  manner  are 
displayed  in  figure  24  for  Lexan*  rolled  homogeneous  armor  steel* 
armor  aluminum,  soft  aluminum  and  lead.  From  static  tests,  we  oan 
fit  our  results  for  all  materials  to: 

v.  '■  *j(f)'6  (3°) 

where  ko  is  different  for  each  material.  A  summary  of  results 
for  a  wide  variety  of  materials  is  shown  in  figure  25. 

The  static  test  results  summarized  in  the  figure  lead 
immediately  to  one  prediction  for  the  dynamic  impact  tests.  One 
of  our  dynamic  tests  consisted  of  dropping  a  .250-inoh  WC  ball  on  a 
target  from  a  height  of  14  ft  to  obtain  a  low  velocity  impact  data 
point.  If  this  impact,  which  we  call  the  "Rafter  Drop  Experiment," 
produced  a  hole  .001-inch  deep  or  greater,  which  was  the  limit  of 
measurement  sensitivity,  corresponding  to  p/d  >  .004  ,  then  we 
concluded  that  the  target  had  been  damaged.  The  Impact  velocity 
corresponding  to  a  14-ft  height  is  28  ft/sec.  Referring  to 
figure  25,  we  draw  a  vertical  line  at  p/d  ■  .004  .  One  would 
predict  that  any  material  whose  V*  1b  greater  than  28  ft/sec  at 

this  p/d  should  show  no  damage,  because  the  energy  could  all  be 

absorbed  elastically.  The  plastic  work  done  at  this  p/d  is 

negligible.  Any  material  with  V#  <  28  ft/sec  at  p/d  ■  ,004  , 

we  would  predict  should  show  surface  damage  at  least  one  mil  deep. 
This  is  what  was  found  experimentally.  The  lead,  soft  aluminum, 
zinc  and  armor  aluminum  showed  surface  damage  in  the  rafter  drop. 
The  steels  and  Lexan  did  not.  The  cast  iron  showed  marginal 
damage  in  the  rafter  drop,  as  would  be  expected  from  its  V#  of 
28  ft/sec  at  p/d  ■  .004  . 
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When  the  static  test  results  are  converted  into  values  for 
E#e  ,  the  curves  of  figure  26  are  obtained.  Lexan,  which  has  a 
low  density  and  high  elastic  energy,  has  a  very  large  E#e  com¬ 
pared  to  the  metals.  E*e  was  found  from  static  tests  to  be  well 
described  for  all  materials  by: 

E.e  ■  k3(l)''75  (31) 

where  kj  is  a  different  constant  for  each  material. 

We  summarize  our  results  in  the  following  way: 


Experimentally,  we  found  V#  -  (p/d) * 6  .  while  theoretically 
we  would  have  expected  roughly:  V#  -  (p/d)*  .  The  additional 

term  (1  -  p/d),7s  will  have  the  effect  of  reducing  the  ex¬ 
ponent  somewhat,  but  the  discrepancy  in  exponents  probably  is  due 
to  our  assumption  that  the  contact  area  does  not  increase  as  force 
is  applied  to  the  ball. 

Similarly,  for  E*e  the  experimental  value  is  (p/d)"’*71  , 
while  theoretically  it  was  predicted  that  E*«  “  (p/d)“*s  .  The 
approximate  agreement  of  the  expressions  and  the  exponents  leadB 
us  to  believe  we  understand  the  mechanism  of  elastic  energy. 

We  find  from  the  static  tests  that  a  value  for  ki  of  0.62 
provides  a  good  fit  to  all  the  data.  Using  the  empirical  (p/d) 
dependence  and  incorporating  the  theoretical  coefficients,  wo 
synthesize  the  following  formulas  for  V#  and  E«e  • 
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Theoretical  (Synthesis) 

v* " 3,3  ^ 

(32) 

E„  -  1.09  •  10‘3  |^(|)'-75 

(33) 

where  B  is  Brinell  hardness  In  Nt/ma  ,  E  is  Young's  modulus 
in  Nt/ma  ,  pt  is  the  target  density  in  kg/m’  ,  mp  is  the 
penetrator  mass  in  kg  ,  E*e  is  given  in  BTU/lbm,  V*  in  ft/sec. 

In  figure  27  the  theoretical  values  of  the  constant  coeffi¬ 
cients  are  compared  to  the  values  found  in  static  tests.  The 
theoretical  expression  shows  good  rough  agreement  for  a  wide  range 
of  materials.  We  conclude  that  if  one  cannot  perform  the  static 
tests  in  advance  to  measure  kg  and  kj  ,  the  expressions  in 
equations  (32)  and  (33)  will  provide  a  good  approximation. 

As  an  illustration  of  the  importance  of  E*e  for  some 
materials,  refer  to  the  data  for  Lexan  in  figure  28.  The  Em 
obtained  in  high  speed  impact  tests  was  found  not  to  be  constant, 
but  to  become  very  large  at  low  p/d  .  The  solid  line,  which  is 
the  predicted  E#  from  static  tests,  and  which  goes  as  ( p/d ) “ • 75 
for  small  p/d  ,  fits  the  data  very  well,  We  conclude  that  for 
Lexan  the  elastic  energy  is  comparable  to  the  plastic  work  for 
p/d  as  large  aB  2.  Thus,  for  highly  elastic  materials,  E#e  can 
dominate  E„p  , 

Theory  for  the  Plastic  Energy  E#p 


The  plastic  deformation  work  performed  on  the  target  is  de¬ 
scribed  by  EMp  .  The  plastic  deformation  may  be  thought  of  as 
analogous  to  a  Brinell  Hardness  Experiment  (Bee  figure  29a).  As 
the  particle  moves  through  the  material  there  will  be  a  pressure 
on  the  front  face  which  is  F/wrz  where  r  is  the  contact 
radius.  The  force  F  arises  because  the  particle  does  work  in 
the  volume  around  it  by  shearing  the  target  material  and  causing 
it  to  flow.  The  work  done  on  each  small  region  of  target  material 
is  Just  the  flow  stress  cp  times  the  distance  it  moves,  so  the 
total  work  per  deformation* volume  is  proportional  to  the  flow 
stress  op  ,  Using  the  well  known  Prandtl  solution  to  the  defor¬ 
mation  flow  field  one  finds: 

•  p  /  WORK  S 

■  B  ■  3,iop  -  nVCLV'ME '  (3*0 

irr 

where  B  is  the  Brinell  hardness,  Since  E#p  io  also  a  measure 
of  the  plastic  shear  work  per  volume  done  by  a  penetrator,  v/e 
take : 


COMPARISON  OF  ELASTIC  COEFFICIENTS 
FROM  THEORY  S  STATIC  TESTS 
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where  pp  is  the  density  of  the  target  material  and  op  ,1s  the 
measured  flow  stress  at  temperature  T  and  strain  rate  e  of  the 
experiment . 

It  has  been  found  experimentally  and  there  is  theoretical 
basis  for  this  finding;  that  the  constant  temperature  flow  stress 


may  be  written: 


o«(T,£5  ■  o. 
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(36) 


This  result  comes  from  the  well-known  "velocity-temperature’  or 
"temperature-strain  rate"  interrelationship.  The  import  of  the 
formula  is  that  a  hardness  test  done  at  a  higher  strain  rate 
(assuming  no  heating  effects)  will  be  equivalent  to  one  done  at 
lower  temperature.  $  Is  typically  .018  and  e  ■  lO'Vaeo  , 
where  T  is  expressed  in  °K  .  0 

At  the  very  high  strain  rates  of  impact  the  heat  dissipated 
in  the  material  by  the  flow  work  cannot  be  conducted  away  fast 
enough  to  allow  one  to  use  ordinary  Brinell  hardness  tests  to 
evaluate  E„  .  A  Brinell  hardness  experiment  is  low  enough  so 
any  heat  generated  has  time  to  diffuse,  so  all  the  material  re¬ 
mains  at  one  temperature;  but,  in  an  impact  test,  the  high  strain 
rates  of  shear  near  the  penetrator  can  cause  large  local  tempera¬ 


ture  rises  and  these  will  in  turn  change  the  flow  stress 
those  regions. 


in 


To  account  for  this  phenomenon  when  using  Eq.  (35)  we  must 
solve  for  the  temperature  rise  caused  by  the  flow  work.  The 
equation  for  this  is 

aPT5f +  +  ‘u  ^  (37! 


where  the  total  stress 
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i0  iB  the  bulk  coefficient  of  thermal  expansion.  This 
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and 

equation  describes  the  temperature  rise  in  a  small  control  volume 
of  target  material  (see  figure  29b)  as  it  passes  around  the  pene¬ 
trator.  Cp  is  the  heat  capacity  of  the  target  material.  The 
first  term  on  the  right  accounts  for  adiabatic  compressional 
heating  of  the  control  volume,  kV*T  accounts  for  any  heat  flow 
into  or  out  of  the  volume,  and  the  last'  term  expresses  the  shear 
heating  of  the  small  volume  of  material. 

Now  we  shall  argue  that  for  our  impact  experiments  the  first 
two  terms  on  the  right-hand  side  may  be  neglected  compared  to  the 
third  term.  The  third  term,  which  is  the  flow  stress  times  shear 
strain  rate  is:  - 

aui  ,f 


V 


whore  a  is  the  radius  of  the  penetrator  and  V  its  velocity. 
The  adiabatic  congressional  term  is  Just  a  T  times  the  change 
of  pressure  with  time  which  is:  p 


,TH-  °[v? 


(op  +  pV4 


For  the  range  of  impact  velocities  of  interest,  rV*  <  o«  . 
In  addition,  ap  -  10-*/°C  and  T  -  10+7°C  so: 


.T  M  ■  0  [10"2°F  s] 


and  is  therefore  negligible  relative  to  the  shear  heating.  And, 
finally, 

*v2t  -  o(*§) 

CL 

For  materials  where  the  strain  rates  are  so  high  chat : 

V  kT 
°F  a  >> 

a 

a  simple  theory  may  be  developed.  For  example,  in  the  case  of 

it*"*  .  v 

copper,  which  has  a  very  large  k  ,  ^  •  10"*  a y  £  for  typical 

a 

Impact  tests.  Equation  (22)  may  then  be  written: 

‘  °F  i 


pcpdT  s  fi°F  r 


where  z  measures  the  distance  traversed  by  a  material  clement, 
as  in  figure  29b. 

We  assume  that  the  temperature  dependence  of  the  flow  stress 
can  be  described  by:  ,m  m  v  /mm  \ 


/T_T\  / T-T0 


Op.(T)  ■  op(To)e  '  1  op  e  '  1  m/  (*U) 

*At  high  velocities,  -  VgouND  in  the  target,  compressions! 
heading  will  become  important,  since  in  this  region  pV2  >>  op 

When  V  -  VS0UND  ,  then  PV2  -  E  ~  102ap  ,  or  apT<$|  “  op  J  . 
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This  is  in  good  agreement  for  all  materials,  where  T;n  is  the 
melting  temperature  and  a  <  <  .4  .  Then  Eq.  (40 )  can* be 

written 


pCpdT 


/T-Tc’\ 


(42) 


This  may  be  integrated  to  give 
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(43) 


Ae  the  volume  of  target  material  progresses  a  distance  z  along 
its  flow  path,  the  local  flow  stress  decreases  according  to  this 
relation. 


The  average  level  of  flc-  stress  in  the  deforming  medium  in 
the  vicinity  of  the  penetrator  can  be  written 


<o,> 

J.‘ 
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(44) 


which  is  the  volume  average  of  op  ,  If  (Z2-Z1)  is  taken 
f2a  ,  where  a  is  the  ball  radius,  then: 
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equal  to 


(45) 


Note  that 


<oF> 


only  really  depends  on  one  parameter: 
1 11  2 


(46) 


and  the  measurable  static  properties  of  the  material,  using 
Eq.  (36)  to  account  for  strain  rate. 


We  reason  that  a  is  roughly  a  constant  for  the  following 
reasons: 

1.  oi  does  not  vary  too  much  (.2  <  a]_  <  .4) 

2.  f^  is  a  measure  of  the  gradient  of  the  shear  flow  field. 
If  6  is  a  characteristic  thickness  of  the  shear  flow 
layer,  then  f]_  -  a/6  . 

3.  f2  is  proportional  to  the  length  of  the  integration 
region.  The  thickness  of  the  shear  layer  is  a  measure 

of  the  length  over  which  the  target  material  is  softening. 
Then  f2  ~  6/a  . 
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Thus,  for  the  purpose  of  constructing  an  approximate  theory,  we 
assume  a  is  a  constant, 

Using  a  which  we  assume  to  be  constant,  we  substitute 
<0p>  of  Eq.  ( 45 )  for  op  of  Eq .  (35),  obtaining: 


E»p 


C5oC  T  In 
2  pm 


(aT^i  * 


c ) 


(47) 


Correlating  this  theoretical  result  with  the  data  from  our 
impact  experiments,  we  arrive  at  values  for  the  constants  a  and 
C2  .  In  this  equation  at  the  current  time  we  use: 

a  ■  0.05 
C2  ■  6.64 

We  have  found  that  this  formula  predicts  E#p  for  a  wide  variety 
of  materials  to  about  s  16%  , 

In  figure  30  the  experimental  value  of  E#c  is  compared  to 
the  theoretical  result  from  Eq.  (47).  The  solid  squares  plot  the 
ratio  of  experimental  to  theoretical  E#p  for  materials  as  dis¬ 
similar  as  Lexan,  steel,  and  lead.  All  the  points  fall  within  the 
error  band  of  ±  1636  .  As  a  means  of  comparison,  the  open  circles 
plot  the  ratio  of  experimental  Ewp  to  the  heat  of  fusion  Hf  of 
the  target  material.  Hf  has  sometimes  been  used  as  an  indicator 
of  E-p  .  One  sees  in  figure  30  that  Hf  works  well  for  Home 
materials  but  not  for  lead,  copper  or  aluminum.  In  addition, 
there  is  no  Hf  for  Lexan.  In  comparison,  our  formula,  Eq.  (47), 
predicts  accurately  the  E#p  for  all  these  material?  as  well  as 
others  presented  in  this  report. 

The  total  E*  of  a  target  material  is  the  sum  of  the  plastic 
part  and  the  elastic  part: 


(48) 

E*P  *a  a  constant  vs.  p/d  ,  while  E#e  -  (p/d)”*  78  .  For 
most  of  the  materials,  particularly  the  pure  metals  presented  in 
Chapter  3,  the  elastic  contribution  is  quite  small,  so  E#t  is 
almost  constant  with  p/d  .  A  typical  example  of  this  is  zinc, 
shown  in  figure  31.  Just  the  reverse  holds  true  for  Lexan, 
figure  28,  where  -E*e  >  E#p  for  p/d  <  2  . 

For  the  entire  range  of  materials  tested,  the  results  of 
which  were  presented  in  Chapter  3,  the  theoretical  predictions  ob¬ 
tained  from  E#p  were  found  to  be  in  good  agreement  vjith  the  data 
over  a  wide  range  of  p/d  . 

The  true  test  of  a  theory  is  whether  it  can  accurately  pre¬ 
dict  the  outcome  of  an  experiment  in  advance  of  the  experiment. 

In  figures  32,  33,  and  34  we  present  predictions  for  three 
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EXPERIMENT  vs.  THEORY  FOR 


Figure  30 


ARMOR  STEEL  TARGET 


Figure  33 


BORON  CARBIDE  TARGET 


Figure  34 
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different  materials  made  several  weekB  before  the  impact  tests 
were  carried  out.  Figure  32  presents  the  prediction  and  experi¬ 
mental  data  for  armor  aluminum.  The  data  clearly  show  the  pre¬ 
dicted  elastic  behavior  at  low  p/d  as  well  as  the  asymptotic 
approach  to  E„p  for  large  p/d  .  The  overall  agreement  between 
experiment  and  theory  is  within  better  than  20%.  In  figure  33  is 
presented  the  data  for  armor  steel  and  the  theoretical  curve,  also 
predicted  in  advance  of  experiment.  Again,  the  elastic  region  is 
discernible,  and  the  overall  agreement  with  theory  is  15X  or  better. 

A  test  of  great  interest  to  us  was  that  of  boron  carbide,  B.'jC  , 
since  our  theory  for  E#  predicts  that  all  the  light,  hard 
ceramics,  particularly  B4C  ,  should  have  very  high  Em's  ,  many 
times  that  of  armor  steel  or  armor  aluminum.  Our  prediction  of 
Egp  for  Bi|C  is  1197  BTU/lbm,  compared  to  around  200  BTU/lbm 
for  steel  and  aluminum  armor.  In  addition,  the  elastic  contri¬ 
bution  is  large  at  low  p/d  .  Using  the  published  value  of  hard¬ 
ness  of  2800  kg/mm*  ,  published  values  of  E  ,  p  ,  T  ,  and  Cp  , 
we  obtained  from  our  formulas  the  prediction  plotted  as  a  solid 
line  in  figure  31**  Our  first  three  data  points,  which  were  ob¬ 
tained  afterward,  are  also  plotted.  Note  the  vertical  axis  is 
logarithmic  with  values  of  10,000  BTU/lbm  obtained  at  low  p/d  , 
and  approaching  1,000  BTU/lbm  at  p/d  -  A  .  Because  the  targets 
we  used  were  -  .4"  thick,  baokfaoe  effects  beoame  important  at 

larger  p/d  ,  quite  possibly  reducing  the  experimental  value  of 

E#  .  Thicker  targets  and  higher  velocity  experiments  are  required 
to  check  our  theory  at  deeper  penetration.  Nevertheless,  it  1b 
clear  that  our  theory  is  rather  successful  in  predicting  the  ex¬ 
tremely  good  armor  properties  of  boron  carbide. 

Thus,  we  believe  we  understand  the  mechanisms  involved  in 
armor  and  we  have  derived  a  theory  which  is  good  at  the  extremes 
of  very  soft  materials  (lead)  to  very  hard  materials  (B4C),  from 
very  elastic  materials  (Lexan),  to  very  inelastic  materials  (zinc). 
For  lightweight  armor  what  is  required  is  a  high  E*  .  Using  the 
formula  of  Eq.  (47)  for  E#p  we  have  compiled  a  table  of  69 

materials  for  which  hardness  data  and  other  parameters  could  be 

obtained,  presented  in  Appendix  I.  From  this  we  r.ummarize  in 
Table  3  those  materials  with  the  highest  E#p ,  which  should  there¬ 
fore  make  the  best  lightweight  armor. 

On  the  list  appear  a  number  of  materials  which,  it  is  known, 
are  good  armor  materials,  such  as  Be2C  ,  B^C  ,  BeO  ,  TiBj  ,  AI0O3 
etc.  Silicon,  hard  chromium  and  other  materials  are  included  at 
the  end  of  the  table  for  comparison. 


The  ranking  presented  here  is  based  on  E«  alone,  which 
assumes  no  deformation  of  the  penetrator.  If  this  in  included, 
the  pressure  developed  at  the  front  face  of  the  penetrator  also 
becomes  a  consideration.  One  wants  the  highest  practical  pressure 
at  the  leading  face  in  order  to  plastically  deform  the  penetrator. 
For  low  velocities  the  pressure  is  Just  PtE*p  -  Thus,  for  a 
given  E#p  ,  one  desires  the  armor  with  the  highest  density.  For 
example,  T1B2  ,  although  it  has  an  E#p  20*  lower  than  BjjC  , 
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has  78S  higher  density  and  may  make  a  better  armor  because  it 
deforms  the  penetrator  more.  A  program  is  currently  underway  to 
include  In  the  above  analysis  the  effects  of  deforming  penetrators. 


PROMISING  ARMOR  MATERIALS 
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CHAPTER  5 


Conclusions 


A  wide  range  of  target  materials  has  been  tested  In  the 
A , R, A , F ,  Impact  Faoility  and  their  Impact  properties  have  been 
evaluated  using  the  Integral  theory  of  impact.  The  results 
demonstrate  that  good  correlation  can  be  achieved  between  theory 
and  data  using  the  two  characteristic  material  properties  E«  and 
V#  (or,  alternatively,  E„p  and  E#e). 

A  theory  has  been  developed  which  relates  these  two  parameters 
to  more  fundamental  material  properties.  The  agreement  between 
thoory  and  the  experimentally  deduced  value  of  E*  and  V*  is 
very  good.  The  theory  predicts  a  plastic  contribution  to  EH 
which  is  roughly  constant  for  a  given  material  and  an  elastic 
contribution  which  is  a  function  of  the  depth  of  penetration.  The 
elastic  contribution  is  small  for  most  metals  but  for  some 
materials  suoh  as  polycarbonate,  it  is  the  dominant  effeot  for 
small  penetrations. 

The  results  demonstrate  that  for  present  day  armors  E#  is 
typically  200  BTU/lbm.  The  theory  suggests  many  candidate  armor 
materials  with  vastly  superior  E„  capability.  One  of  these 
materials,  boron  carbide,  has  been  tested  and  a  value  of 
Em-  1200  BTU/lbm  has  been  obtained  -  a  value  in  good  agreement 
with  the  theory, 

It  is  now  possible  to  utilize  the  knowledge  obtained  in  this 
program  to  optimize  the  design  of  armor  (or  penetrator)  systems. 
This  is  the  object  of  the  next  portion  of  the  program. 
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